Magnetization plateaus as insulator-superfluid transitions in quantum spin systems 
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We study the magnetization process in two-dimensional S = 1/2 spin systems, to discuss the 
appearance of a plateau structure. The following three cases are considered: 1) the Heisenberg 
antiferromagnet and multiple-spin exchange model on the triangular lattice, 2) Shastry-Sutherland 
type lattice, (which is a possible model for SrCu2(B03)2,) 3) 1/5-depleted lattice (for CaV409). We 
find in these systems that magnetization plateaus can appear owing to a transition from superfluid 
to a Mott insulator of magnetic excitations. The plateau states have CDW order of the excitations. 
The magnetizations of the plateaus depend on components of the magnetic excitations, range of the 
repulsive interaction, and the geometry of the lattice. 

PACS numbers: 75.60.Ej, 75.10.Jm 



In some one-dimensional spin systems, spin-density- 
wave states with finite spin gap appear under a finite 
magnetic field accompanying plateau structures in the 
magnetization process. Magnetization plateaus were ob- 
served in some quasi one-dimensional materials.! The- 
oretical arguments clarify that the appearance of the 
plateau is explained by an insulator-conductor transition 
of magnetic excitationsfl In two- or higher-dimensional 
systems, magnetization plateaus have been also ipund in 
both theoreticaflij- and experimental studiesfliS In this 
paper, we propose a rather general picture that these two- 
dimensional plateaus are formed owing to field-induced 
insulator-superfluid transitions of magnetic excitations. 
To demonstrate how it works, we discuss three examples 
in details. 

The first example is a family of antiferromagnets on 
a triangular lattice. For the S = 1/2 antiferromagm 
on a triangular lattice (AFT), Nishimori and Miyashit. 
found a magnetization plateau at ?7i/msat = 1/3, which 
comes from the appearance of a collinear state with three 
sublattices, i.e., the so-called "uud" state. This plateau 
was actually observed in AFT materials like C6Eu(Ref. |g) 
and CsCuCl3(Ref. M). Recently in a multipla-spin ex- 
change (MSE) model, which is a possible modelH for solid 
•^He films, a magnetization plateau was predictedEI at 
m/msat — 1/2. In this case, the plateau is attributed 
to the formation of a similar collinear state but with 
four sublattices. The magnetization processes of these 
systems have been studied extensively and here we just 
attempt interpreting the known results to test the new 
picture. 

We take as the second example the S = 1/2 Heisenberg 
antiferromagnet (HAF) on the Shastry-Sutherland lat=. 
tice (Shastry-Sutherland model, hereafter. See Fig. 0),t2l 
which is known to have-jan exact dimer ground state. Re- 
cently Kageyama et al.El found that SrCu2(B03)2 realizes 
a lattice structure equivalent to that discussed in Ref. |l^ 
and that it has a gapful ground state. The magnetization 
measurements show plateaus at m/msat = 1/8 and 1/4. 



The last is the S = 1/2 HAF on the 1/5-depleted square 
lattice (Fig. g), which includes a model Hamiltonian for 
CaV409. In this system, the. plaquette singlet state is 
realized in the ground statcllil 

In our picture, the plateau states can be regarded as 
Mott insulators of effective magnetic particles; repulsive 
interactions induce various kinds of charge-density-wave 
(CDW) long-range order leading to a finite energy gap in 
particle-hole excitations. Except for the plateau phases, 
magnetic particles are conducting to form supersolid, in 
which superfluidity and CDW coexist, and magnetiza- 
tion increases smoothly. Of course, the charge density 
is translated into the spin (S^) density and superfluidity 
here means long-range order in the direction perpendicu- 
lar to the field. Although this essential picture is common 
to the three examples, the concrete forms of the magnetic 
particles are different. 

In the first example, i.e. AFT and an MSE model, a 
single fiipped spin itself works as the magnetic particle, 
while the triplets on the dimer (J) bonds are the rele- 
vant particles of the Shastry-Sutherland model. We find 
plateaus at m/msat = 1/2 and 1/3. In the third one, i.e. 
the 1/5-depleted square lattice, plaquette triplets behave 
as particles. We predict plateaus at m/rrisnt = 1/8, 1/4, 

1/2. 

The magnetizations where a plateau appears depend 
on (i) the form of the magnetic excitations, (ii) range 
of the repulsion between them, and (iii) the geometry of 
the lattice. Finally we summarize common features on 
properties of the phase transition. 

Spin — magnetic particle: The magnetization plateau 
for the AFTrSystem was found in spin 1/2 anisotropicnl 
and isotropicQ Heisenberg models. The Hamiltonian is 
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7jS!S^)-BJ2S!, (1) 



where the summation runs over all nearest-neighbor pairs 
and B denotes the magnetic field. For 77 > 1, the mag- 
netization curve has a plateau at m/rrisat — 1/3. The 
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FIG. 1. Shastry-Sutherland lattice 
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FIG. 2. 1/5-depleted square lattice. 



ground state in the plateau phase is of a coUinear struc- 
ture with three sublattices, where two of three spins 
direct upward and the other downward. In the other 
phases, magnetic states have non-coUinear structures 
with off-diagonal long-range order (ODLRO). If we in- 
troduce a particle picture, i.e., recognize spin dynami^ 
as induced by the motion of a certain kind of particles ,tHl 
the appearance of plateau is easily understood from sim- 
ple consideration about compressibility of the particles. 
Regarding an up ^in as a hard-core boson and a down 
one as a vacancy,t3 we can rewrite the Hamiltonian (U^) 
as 

H = - ^{{h^^o + h.c) + 277n,nj} - (B + 3J77) ^n„ 



{»,j> 



(2) 



where 6] denotes the creation operator of the hard-core 
boson on site i, and rii the number operator. Since parti- 
cles carry magnetic moment unity, the chemical potential 
H = B + 3Jri is controlled by the magnetic field and the 
/i-dependence of the particle density n corresponds to the 
magnetization curve of the original spin system. 

The hopping term comes from the spin exchange (XY) 
term and the repulsive interaction from the diagonal 
(Ising) part. The anisotropic case, 77 > 1, is mapped to 
the strong coupling (i.e. strong repulsion) region of the 
corresponding boson system. The particle-hole transfor- 
mation converts the system into that of holes with repul- 
sion of the same strength; in the strong coupling limit, 
the ground state at the filling n = 2/3 (m/msat — 1/3) 




FIG. 3. Spin-density-wave order in AFT (a) and MSE 
model (b). Black (white) circles denote down (up) spins. 



has the density wave long-range order, with the three- 
sublattice structure shown in Fig. ||(a). Due to the 
repulsive interaction, this state is incompressible, i.e. 
dn/dii = 0, and density-fluctuation energy has a finite 
gap above the ground state. Except for the filling(s) 
n = 2/3 (and 1/3), there are vacancies and hence parti- 
cles are mobile (conducting) . Since the particles obey bo- 
son statistics and the system is uniform, the system pre- 
sumably shows superfluidity. There is perfect correspon- 
dence between the above consideration and the previous 
results;0u the insulating CDW state with dn/dfj. = 0, is 
consistent with the spin collinear state, where suscepti- 
bility is vanishing. On the other hand, superfluidity of 
bosons corresponds to non-coUinear ODLRO of spins. 

The particle density where CDW stabilizes depends 
on the range of repulsive interactions. To see this, we 
next discuss MSE model with four-spin exchange on the 
triangular lattice, where repulsion acts further than in 
the Heisenberg model. The Hamiltonian is given by 



n = jJ2, 



Kj2hp~Bj2'^t 



(3) 



;j> 



where (Ti denote Pauli matrices. The second summa- 
tion runs over all minimum diamond clusters and hp 
is the four-spin exchange hp — 4(P4 + P4^^) — 1 with 
the ring permutation of four spins P4. It was shown 
that three- and four-spin exchange interactions are very 
strong in two-dirpensional solid "^He due to strong quan- 
tum fluctuations.Ll Theoretically a magnetization plateau 
was founcB at m/msat — 1/2 instead of m/msa,t = 1/3. 
In the particle picture, bosons feel the following two-body 
repulsion 



V = 4( J + 5K) ^ Tiirij +4:K ^ n.^rij . 



(4) 



ihj) 



(i.i) 

ejv.iv.jv. 



The repulsive interaction acts in both nearest- and next- 
nearest-neighbor sites. Figure H(b) shows the region 
where the interaction works. Because of the range of re- 
pulsion, the particles can solidify at the density n — 1/4. 
(Note that the solidification occurs only if the repulsion 
overcomes the effect of the hopping term.) This insulat- 
ing phase at n = 1/4 corresponds to the magnetization 
plateau at m/msat = 1/2 in the original spin system. 
The previous numerical result in Ref. H on the ground 
state of the plateau phase is consistent with the CDW 
order shown in Fig. @(b). 
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FIG. 4. Two-body repulsive interactions up to the 3rd or- 
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Dimer triplet: When specific pairs of two spins are 
more strongly coupled than to the others by an anti- 
ferromagnetic interaction, the dimer singlet state real- 
izes in the ground state. Under a weak magnetic field, 
S^ = 1 triplets on the dimer bonds are dominant exci- 
tations. Several types of repulsive interactions between 
these dimer triplet excitations can induce various insu- 
lating phases and thereby yield magnetization plateaus. 

Forj-cpncreteness, we discuss the Shastry-Sutherland. 
modefij shown in Fig. nl The exact dimer ground stateE2l 
realizjfis for J' /J < 0.69. (Ref. [l^.) Recently Kageyama 
et al.Q found this lattice structure in SrCu2 (603)2 and 
observed the magnetization plateaus at m/msat = 1/8 
and 1/4. Because of the special struoture of the lattice, 
a triplet excitation is almost localizedll3 Considering the 
dimer triplet state with S'^ = 1 as a particle (a hard-core 
boson by definition) and the dimer singlet as a vacancy, 
we derive an effective Hamiltonian for it using the pertur- 
bational expansion from the J' = limit. The expansion 
is performed up to the 3rd order in J' / J from degener- 
ate states with a constant number of dimer triplets. The 
effective Hamiltonian up to the 2nd order is 



J'- 



J 



^-U-^-V E-+ T + ^ E 



J' J'- 
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niUj 



(^i) 



J' 



+ jy ^{[bl{bt+el ~ &j-el) + h.c.]{ni^e2 - ni+e2) 



4J 



ieA 



-f2ni+e2(l - ni)n.i^e2 + (&l+el^*-el + h.c.)ni} 

+ T7E{ ^l^^2 }, 



(5) 



ieB 



where i(j) runs over an effective square lattice consisting 
of dimer bonds (both horizontal and vertical) and hor- 
izontal (vertical) ones belong to A (B) sublattice. The 
full form of the effective Haroiltonian up to the 3rd order 
will be published elsewhere.til The derived Hamiltonian 
does not have the one-particle hopping term (as was al- 
ready reported in ref. |l^), but contains many correlated- 
hopping processes, where an effective hopping of a parti- 
cle is mediated by another one. This is one of our main 
observations. Most 3rd-order terms concern the corre- 
lated hopping. Longer-range repulsions between parti- 
cles appear from higher-order perturbations. Diagonal 




FIG. 5. Magnetization process of the Shastry-Sutherland 
model with J' /J = 0.45. 



repulsive interactions up to the 3rd order in J' /J are 
shown graphically in Fig. 0. The resulting Hamiltonian 
does not have 90° rotational invariance, since the lattice 
structure has low symmetry, and this may lead to highly 
anisotropic CDW states. 

We study the effective Hamiltonian in the classical 
limit. To this end, we map the hard-core boson system 
to the 5* = 1/2 quantum spin system and then approxi- 
mate the spin-1/2 by a classical unit vector. We search 
for the ground state with large sublattice structures (e.g. 
a stripe-like one with 6-sublattice) both with the mean- 
field approximation and a Monte Carlo method by de- 
creasing temperatures gradually. The evaluated magne- 
tization process is shown in Fig. ||. Note a clear dif- 
ference between the high- and low-field region. There 
appear plateau structures at m/msat ~ 1/2 and 1/3. 
The plateau states have CDW long-range orders shown in 
Fig. |6[ Configurations reahzed for m/rrisat = 1/2 and 1/3 
correspond to perfect closed packings provided that par- 
ticles avoid repulsion from 1st- and 2nd-order perturba- 
tion, respectively. The plateau at m/msat = 1/2 appears 
only in the region < J' / J < 0.50 and, for large J' /J the 
CDW is destroyed by the correlated hoppings, which are 
dominant in the higher order terms. The correlated hop- 
pings are so efficient also at large particle density that 
any plateau does not appear for 1/2 < m/nisat < 1. 
Below m/msat — 1/3, the correlated hoppings occur 
rarely, because of a low particle density. The observed 
1/4-plateau (and 1/8-plateau) of SrCu2 (603)2 may be 
formed by weak longer-range repulsions which are not 
taken into account in the present study. Recently Miya- 
hara and Ueda discussed semi-phcnomcnologically that 
the 1/4-plateau ptaight come from a CDW state with a 
stripe structure.llj In our approach, the repulsive inter- 
action relevant to the stripe CDW may come from the 
higher-order terms in perturbation, otherwise from other 
spin interactions that are not considered in the Shastry- 
Sutherland model. This remains to be a future problem. 
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FIG. 6. Spin configurations at magnetization plateaus at 
m/rrisat = 1/3 (a), 1/2 (b). Black bonds denote dimer triplet 
excitations. 



Plaquette triplet: When the interactions or a special 
geometry of the lattice allows four-spin plaquettes, in 
each of which four spins are coupled more strongly than 
to the others, individual plaquettes form singlets in the 
ground state. The triplet states with 5'^ = 1 on plaque- 
ttes are dominant excitations in a weak magnetic field. 
The insulator-conductor transition of these excitations 
can take place thereby producing magnetization plateaus. 

An example of the plaquette singlet ground states is 
seen in— the 5 = 1/2 HAF on the 1/5-depleted square 
lattice,liil which includes a possible model for CaV409 as 
a special case. The lattice is shown in Fig. |[ In the 
isolated plaquette limit Ji — J2 = 0, a, trivial plateau 
already appears at m/nisat — 1/2 for J < B < 2J 
(see Fig. 0), where every ^plaquette is in the triplet ex- 
cited state with S^ = 1.113 When m/msat < 1/2, the 
triplet excitations (particle) tend to hop if Ji and J2 
are turned on, and at specific (commensurate) values of 
m/jTisat they can show insulator-conductor transitions as 
a consequence of the competition between the hopping 
and the repulsive interaction. Above m/msat — 1/2, 
a plaquette quintuplet {S — 2) with S^ = 2 behaves 
as a particle and can show magnetization plateaus be- 
tween 1/2 < m/msnt < 1- In the following, we focus 
on a weak magnetic-field region, which corresponds to 
the magnetization < m,/m,sa.t < 1/2. Regarding the 
plaquette triplet excitation with S*^ = 1 as a particle 
and the singlet state as a vacancy, we derive the effective 
Hamiltonian of the particle by the 2nd-order perturba- 
tion around the limitJi — J2 = 0. (The explicit form will 
be shown elsewhere.Lj) The repulsive interactions range 
from a plaquette to its nearest- and next-nearest neigh- 
bors. If parameters satisfy Ji ~ 2J2, the hopping term 
is weak and hence the system is in the strong coupling 
regime. Then we may expect that the triplets crystallizes 
and the magnetization plateaus appear at m,/m,sa.t — 1/8 
and 1/4. The mean-field approximation of the effec- 
tive Hamiltonian indeed shows magnetization plateaus at 
"m/msat = 1/8 and 1/4 (see Fig. R). They come from the 
insulating phases with CDW long-range order of a square 
structure. A remark is in order here about the relevance 
of our results to CaV40g. Quite recently, it was shownEZi 
that plaquettes of another type (metaplaquettes) consist- 
ing of J2-bonds play the main role in CaV409 contrary 




FIG. 7. Magnetization process up to m/msat ~ 0.5 in the 
1/5-depleted square lattice with Ji/J = 1 and J2/J = 0.5. 
The dotted line shows the case Ji = J2 ~ 0. 

to earlier studiesEJ. The physics is, however, almost the 
same also in this case; the metaplaquette excitations be- 
have like particles and, J- and Ji-bonds induce hopping, 
and so on. Xhe detailed results will be reported in a 
longer paper £f 

Common features: To conclude this paper, we discuss 
a few features shared by the three examples. Accord- 
ing to an analogy to many-particle theories, a plateau 
state corresponds to a CDW insulating state and gapless 
ones to supersolids. As the plateau state collapses by in- 
creasing the applied field, superfluidity appears, whereas 
CDW exists in both phases. Let us consider the case of 
2nd order transition, where the magnetization changes 
continuously. Assuming that the on-set of superfluidity 
is well described by the effective Hamiltonian of the ordi- 
nary bosons with a short-range repulsion for low energies, 
we conclude this transition is of the dynamical exponent 
z = 2;E3 magnetization increases linearly like \H — Hc\ 
apart from possible logarithmic corrections. (Note that 
the form is quite different from that in ID.) 

Note added in proof: K. Onizuka et at. recently 
observed a clear 1/3 plateau in SrCu2(B03)2, which we 
had predicted in this paper and had argued to be of a 
stripe structure. 
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